Typeset with jpsj3.cls <ver.l.O> 



Full Paper 



o 

(N 



(N 



•4— > 

CO 



i 

c 

o 
o 



> 
00 

m 
oo 
m 



OS 

o 



Quantum Phase Transitions in Alternating-Bond Mixed Diamond Chains with Spins 
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We investigate the mixed diamond chain composed of spins 1 and 1/2 when the exchange 
interaction is alternatingly distorted. Depending on the strengths of frustration and distortion, 
this system has various ground states. Each ground state consists of an array of spin clusters 
separated by singlet dimers by virtue of an infinite number of local conservation laws. We deter- 
mine the ground-state phase diagram by numerically analyzing each spin cluster. In particular, 
for strong distortions, we find an infinite series of quantum phase transitions using the cluster 
expansion method and conformal field theory. This leads to an infinite series of steps in the 
behavior of Curie constant and residual entropy. 

KEYWORDS: mixed diamond chain, bond alternation, spin cluster, infinite series of transitions, quantum 
phase transition 



1. Introduction 

Frustration plays a crucial role in low-dimensional 
quantum magnetism. 1 ' 2 ) It not only drives magnetically 
ordered states into disordered states by enhancing quan- 
tum fluctuation but also induces magnetic moments from 
the disordered phase. Remarkably, there exist a class of 
models whose ground states are exactly written down as 
spin cluster states (SCSs) because of frustration. A SCS 
is a tensor product of exact local eigenstates of cluster 
spins; a dimer state is a special case. 

One of the well-known examples is the Majumdar- 
Ghosh model, which is a spin 1/2 antiferromagnetic 
Heisenberg chain with next-nearest-neighbor interaction 
whose magnitude is half of the nearest-neighbor inter- 
action. 3 ' This ground state is a prototype of sponta- 
neously dimerized phases in one-dimensional frustrated 
magnets. Many corresponding materials are also found 
as listed in ref. 4. Another well-known example is the 
Shastry- Sutherland model 5 ) for which a corresponding 
material was synthesized 18 years after its theoretical 
prediction. 6 ' 7 ' The ground states of these models are, 
however, nonmagnetic. Namely, the possible magnetic 
(quasi-) long-range order in the unfrustrated counter- 
part of these models is destroyed by the enhancement 
of quantum fluctuation by frustration. 

One of the authors and coworkers investigated a dia- 
mond chain consisting of the same kind of spin, i.e., a 
pure diamond chain (PDC), as an exactly treatable frus- 
trated model. 8 ' They found that, in the spin-1/2 case, 
it has two different SCSs as the ground state depending 
on the strength of frustration. One of them is the non- 
magnetic phase accompanied by the spontaneous trans- 
lational symmetry breakdown (STSB) and the other is a 
paramagnetic phase without STSB. It is also found that 
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this model has a ferrimagnetic ground state in a less frus- 
trated region. 

Modifications of the PDC have been examined by 
many authors. Among them, the PDC with distortion 
has been thoroughly investigated by numerical meth- 
ods. 9 ~ n ) It is found that a natural mineral azurite con- 
sists of distorted PDCs with spin 1/2. The magnetic 
properties of this material have been experimentally 
studied in detail. 12 ' 13 ' Other materials have also been 
reported. 14,15 ' The effects of the 4-spin cyclic interac- 
tion have recently been investigated by Ivanov et aZ. 16 ' 
It should be remarked that the diamond chain structure 
is one of the simplest structures compatible with the 4- 
spin cyclic interaction. The finite temperature properties 
of diamond chains consisting of Heisenberg bonds and 
Ising bonds are investigated exactly. 17 ' 18 ' The thermo- 
dynamic properties of a similar classical model with a 
hierarchical structure are also studied. 19 ' The behavior 
of the residual entropy of this model is similar to that of 
the diamond chain. 

In our previous work, 20 ' we introduced another ver- 
sion of the diamond chain that consists of two kinds of 
spins, and named it the mixed diamond chain (MDC). 
Although the MDC is a simple extension of the PDC, it 
is a model belonging to a class different from that of the 
PDC in that its ground state is the Haldane state in the 
weak frustration region. We have investigated the MDC 
consisting of spins 1 and 1/2 depicted in Fig. 1(a) in spe- 
cial detail. 20 ' 21 ' Since the MDC has an infinite number 
of local conservation laws like the PDC, a typical ground 
state is a SCS consisting of spin clusters each carrying 
spin 1. A series of quantum phase transitions take place 
between five ground-state phases with different periodic- 
ities and with or without a STSB. Except in the Haldane 
phase, the MDC has macroscopically degenerate ground 
states, because each cluster has three degenerate ground 
states with spin 1. The SCS structures of the ground 
states arc reflected in the characteristic thermal proper- 
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Fig. 1. Structure of (a) the mixed diamond chain and related 
models: (b) the dimer plaquctte chain, (c) the frustrated ladder. 
The large open circles are spin-1 sites, The small open circles 
and filled circles are spin-1/2 sites. 

ties, as reported in rcf. 21. 

The MDC is related to other important models of 
frustrated quantum magnetism. The dimcr-plaquette 
model 22 ' 23 ) shown in Fig. 1(b), which may be regarded as 
a one-dimensional counterpart of the Shastry-Sutherland 
model, reduces to the MDC if the horizontal bond in Fig. 
1(b) is strongly ferromagnetic. The ground state of the 
dimer-plaquette model in this region has not yet been 
studied. In the MDC limit, however, it becomes clear 
that the ground state is the Haldane state. In addition, it 
turns out that all eigenstates can be expressed in terms of 
the eigenstates of finite-length spin 1 Heisenberg chains 
for an arbitrary strength of the plaquette bond in the 
MDC limit. 20 ) ' 

The MDC is also related to rung-altcrnating frustrated 
Heisenberg ladders (Fig. 1(c)). If the spin 1 site in MDC 
is decomposed into the sum of two spin 1/2's, the MDC is 
equivalent to the low-energy sector of the spin 1/2 rung 
alternating frustrated ladder with a leg coupling equal 
to the diagonal coupling and one of the rung couplings 
is weakly antiferromagnetic or ferromagnetic. In the ab- 
sence of rung alternation, it is known that this model has 
only Haldane and rung-dimer ground states. 24 ) However, 
in the MDC limit, the intermediate phases appear as ex- 
plained above. 

Thus far, no materials described by the MDC model 
have been found. Nevertheless, synthesizing MDC mate- 
rials is not an unrealistic expectation in view of the suc- 
cess of the synthesis of many low-dimensional bimetal- 
lic magnetic compounds 25 ) and organic magnetic com- 
pounds. In the latter, spin 1 units included in a MDC 
material can be formed as ferromagnetic dimers. 26 ) Gen- 
erally, materials with a lower symmetry have a higher 
chance of being found or synthesized than those with 
a higher symmetry. Therefore, to raise the possibility of 
realizing MDC materials, theoretical predictions on mod- 
ified versions of the MDC model with a low symmetry 
are required. 

The distortion is a realistic modification that lowers 
the symmetry of the MDC. To classify distortion pat- 
terns, we consider the normal modes of a diamond unit. 




(a) (b) (c) 




Fig. 2. Displacement modes of a diamond unit. 

It has 8 degrees of freedom of displacement within the 
diamond plane. Excluding two translations and one rigid 
body rotation, we have five normal modes as depicted in 
Fig. 2. We may also expect that the distorted MDC can 
be realized as a result of the collective softening of this 
normal mode. Two of them ((a) and (b)) break the lo- 
cal conservation laws, which hold in undistorted MDCs. 
Among other distortion patterns ((c), (d) and (e)), (d) 
and (e) do not change the geometry of the original undis- 
torted MDC. The distortion pattern (c) brings about the 
bond alternation in the undistorted MDC to break the 
original reflection symmetry. This type of distortion not 
only raises the possibility of the experimental realization 
of the MDC, but also is of theoretical importance, be- 
cause in this case the eigenstates are exactly expressed as 
SCSs by virtue of the same local conservation laws as in 
the undistorted case. We hence restrict ourselves to this 
case in the present paper. We will show that the bond al- 
ternation distortion does not destroy the exotic phases of 
the undistorted MDC but even makes the ground-state 
phase diagram more exotic than that of the undistorted 
MDC. Since cases (a) and (b) contain different physics 
owing to effective interactions between spin clusters, they 
will be reported in a separate paper. 

This paper is organized as follows. In §2, the model 
Hamiltonian is presented. In §3, the classical limit is dis- 
cussed. In §4, the structure of the ground state is ex- 
plained and the ground-state phases are numerically de- 
termined. In §5, an infinite series of phase transitions 
that occur for strong bond alternation are examined us- 
ing numerical, cluster expansion, and conformal field the- 
ory methods. The last section is devoted to summary and 
discussion. 

2. Hamiltonian 

The alternating bond MDC is described by the Hamil- 
tonian 

H =^{(1 + 6)3^+^) 
i=i 

+ (l-j)( T W +T }^ J+1 + Ari 1 M 9) }. W 

where Si and (a = 1,2) are spin 1 and spin 1/2 
operators, respectively. The number of the unit cells is 
denoted by N. The Hamiltonian includes three types of 
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Fig. 3. Structure of alternating bond mixed diamond chain with 
5 = 1 and r' 1 ' = t' 2 ) = 1/2. A, 1 + 5, and 1 — <5 are the exchange 
parameters, where A and <5 control the strengths of frustration 
and alternating distortion, respectively. 



exchange parameters, namely, 1 + S, 1 — 5, and A; <5 rep- 
resents the strength of the alternating-bond distortion, 
and A controls the frustration as depicted in Fig. 3. In 
the case of 5 = 0, eq. (1) reduces to the Hamiltonian of 
the MDC without distortion. 20 ' 

3. Classical Ground States 

Before analyzing the full quantum system (4), we ex- 
amine its classical version in which the spin operators 
Si and T\ are replaced by classical vectors with the 
length S and r, respectively. We also impose the condi- 
tion r = S/2. By introducing S t = (l + 5)Si + (l-6)S l+1 



and Ti = t 



(i) 



.(2) 



the classical Hamiltonian is ex- 



pressed in the following two forms: 
HCl = \ 2 [t 2Tl + SO 2 - 2 ( 2 - X ) T f - S 2 - AS 2 



2S 2 



A 2 5 2 



(2) 
(3) 



These equations have the same form as those in the ab- 
sence of distortion, 20 ' since the distortion parameter is 
absorbed in Si. 

For A < 2, the expression (2) shows that H cl is mini- 
mized if \Si\ = 25, \Ti\ = S and \T t + ±S t \ = 0. Then 

all the Si's (t[ q ^'s) in the chain are aligned parallel (an- 
tiparallel) to a fixed axis, and the ground state is an- 
tiferromagnetic. This ground state is elastic, since any 
local modification of the spin configuration increases the 
energy. For A > 2, the expression (3) reveals that W cl is 
minimized if \S t \ = 25" and \T t + S;/A| = 0. Hence, S t 



and Si+i are parallel, and rj 1 ^ and t\*' form a triangle 



.(2) 



with Si/X. t} 1 ^ and t\ z> may be rotated about the axis 
of Si and S,+i without raising the energy. Then all the 
S^s in the chain are aligned parallel to a fixed axis, and 
the arbitrariness of the local rotation of rj 1 ^ and t[ 2 ^ is 
not obstructed. Thus, the ground state is ferrimagnetic 
with magnetization (1 — 2/X)SN. 

Thus, we have two classical phases separated by the 
phase boundary A = 2 independent of 5. In fact, 5 is em- 
bedded in Si and does not explicitly appear in the en- 
ergy expressions cqs. (2) and (3). The distortion affects 
the ground state only in the presence of quantum effect. 
As was pointed out, for A > 2, the classical ground-state 
configuration can be locally modified with no energy in- 



.(2) 



crease. This classical situation corresponds to a quan- 
tum situation in which there are an infinite number of 
low-energy states that are made from each other by lo- 
cal modification. Since such quasi-degenerate low-energy 
states may enhance quantum fluctuations, we expect the 
appearance of exotic quantum states for A > 2 in the 
quantum system (4). 

4. Ground-State Phase Diagram 

The Hamiltonian (1) has a series of conservation laws. 
To see it, we rewrite eq. (1) in the form, 



N 



(1 + 6)SiTi + (1-S)TiSi +1 + ±(t1-'± 



(4) 



where the composite spin operators T; are defined as 

Ti^t^+tP (l = l,2,..-N). (5) 



Then it is evident that 



[T 2 ,-H]=0 {l = l,2,---N). 



(6) 



Thus, we have N conserved quantities T x for all I, even 
if we have introduced an alternating-bond distortion. By 
defining the magnitude 2} of the composite spin Ti by 
Tf = Ti(Ti + 1), we have a set of good quantum numbers 
{Ti; 1 = 1,2, N}. Each T x takes a value of or 1. The 
total Hilbert space of eq. (4) consists of separated sub- 
spaces, each of which is specified by a definite set of {!)}, 
i.e., a sequence of and 1. A pair of spins with 7] = is 
a singlet dimcr. A cluster including n successive T; = 1 
pairs bounded by two Ti — pairs is called a cluster-n. 

Thus, the eigenstates of a cluster-n are expressed in 
terms of the spin states of n T;'s with TJ = 1 and 
n + 1 Si's of Si = 1. Hence, a cluster-n is equivalent to 
an alternating bond antiferromagnctic Hciscnberg chain 
(AAFH1) consisting of 2n + 1 effective spins with spin 
magnitude 1 as in the undistorted case. 20 ' 21 ' We follow 
the terminology of refs. 20 and 21 to call a ground state 
consisting of a uniform array of cluster-n's as the dimcr- 
cluster-n (DCn) phase. The phase boundary between the 
DC(n — 1) and DCn phases is given by 

A c (n - 1, n; 5) ={n + l)E G {2n -1,5)- nE G (2n + 1, 5), 

(7) 

where Eq(L, 5) is the ground-state energy of the AAFH1 
with a bond alternation 5 and a length L. These phase 
transitions are of the first order, since they take place as 
level crossings between two eigenstates of the Hamilto- 
nian (1) characterized by different sets of quantum num- 
bers {!)}. 

The width AA(n; 5) of the DCn phase for each 5 is 
given by 

AA(n; 5) = A c (n - 1, n; 6) - A c (n, n + l;5) (8) 
= (n + 1) [^ G (2n -1,5)- 2E G {2n + 1, 5) 



+E G (2n + 3,5) 



(9) 



Equation (9) shows that AA(n; 5) is proportional to the 
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Fig. 4. Size dependence of the ground-state boundary energy 
AE G (L,8) of AAFH1 for (a) < 8 < 8 C and (b) 8 > 8 C ob- 
tained by numerical diagonalization, and that for (c) <5 = 0.2 
and 0.3 obtained by DMRG. 



second difference of E G (2n + 1,<5) with respect to n. 
Therefore, the DCn phase appears with a finite width 
as long as E G (2n + 1,5) is a convex function of n. If 
AA(n; S) is positive for all n, an infinite series of phase 
transitions take place before reaching the DCoo phase. 

We have calculated E G (L, S) by numerical diagonaliza- 
tion for 3 < L < 19 and by the finite-size density matrix 
rcnormalization group (DMRG) for larger L. To extract 
the boundary contribution, we subtract the correspond- 
ing bulk energy from E G (L, 8) as 



AE G (L, S) = Eq(L, 5)-(L- l)e G (oo, 5), 



(10) 



where £g(oo, 5) is the ground-state energy of the infinite- 
size AAFH1 per spin. This subtraction does not influ- 
ence the second difference. The system size dependence of 
AE G (L, 6) is shown in Figs. 4(a) and 4(b) for 3 < L < 19. 
Figure 4(c) shows the data for L < 31 for 5 = 0.3 and 0.2. 
The energy e G (oo,S) is estimated using the infinite-size 



DMRG by measuring the bond energy of the two bonds 
closest to the center of the chain. 

Figure 4 shows that the ground-state energy of the 
AAFH1 is a convex function of L for large S within nu- 
merical data. Therefore, DCn phases are realized for all n 
depending on A. In contrast, for small S, no DCn phases 
with large n appear, since the ground-state energy of the 
AAFH1 is a concave function of L for large L as shown 
in Fig. 4. Actually, it is known that no DCn phases with 
n > 4 appear for S = O. 20 -* In such cases, the direct tran- 
sition from the DCn phase to the DC oo phase takes 
place at 

A c (n, oo; 5) = E G (2n + 1,6)- 2(n + l)e G (oo, 5), (11) 

if A c (n, oo; 5) < A c (n, n + 1; 5). 

The phase diagram is shown in Fig. 5 using these val- 
ues. At 5 = S c ~ 0.2598, where the Haldane-dimer phase 
transition takes place in the infinite-size AAFH1, 27 ) the 
convergence of the infinite-size DMRG becomes worse. 
Hence, we employed extrapolation from the finite-size 
DMRG results to determine e G (oo,6). For S < 5 C , the 
DCoo ground state corresponds to the uniform Haldane 
phase, while it corresponds to the dimer phase for S > S c . 

The phase boundary between the DC0 and DC1 
phases is analytically obtained. Obviously, E G (1,8) = 0, 
and E G (3,6) satisfies the following eigenvalue equation 
within the subspace with total spin 1: 

Eg{3, S) 3 + 4E G (3, 6f + (3 - 7<5 2 )^ G (3, 5) - 16<5 2 = 0. 

(12) 

Therefore, A c (0, 1) satisfies 

A c (0, 1; 5f - 4A C (0, 1; Sf + (3 - 7<5 2 )A C (0, 1; 8) 

+ 16<5 2 =0. (13) 

This relation is plotted in Fig. 5 by the thick dotted line. 
For 5 ~ 1, eq. (13) implies A c (0, 1; 5) ~ 4 - 2(1 - 5). 

5. Infinite Series of Quantum Phase Transitions 

5.1 Almost decoupled limit: 5 ~ 1 

At 5 = 1, an AAFH1 with a length of 2n + 1 
is decoupled into n independent dimers and a single 
spin 1. Therefore, we can employ the cluster expansion 
method 28 ^ to estimate E G {L,5). The lowest-order non- 
vanishing contribution to AA(n; 8) is 0((1 — S) 2n ) and is 
written as 

A\(n;d)~A n (l + 5) ° 



l + S 



~ ~ n (i _ 5) 



In 



(14) 



The factor A n is given by 



A 1 


= 4/3, 


A 2 


= 2.305555556, 


A 3 


= 4.501164348, 


A 4 


= 8.794932620, 


A 5 


= 17.01798127, 


A 6 


= 32.60950955, 


A 7 


= 61.96667474, 
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Fig. 5. Ground-state phase diagram. Small open circles are the 
critical points A c (n, n — 1; <S) and large open circles are the crit- 
ical points A c (oo,n;<5). The solid lines are guides for the eye. 
For n > 4, phase boundaries are only shown by solid lines to 
avoid complication. The values of A c (0, 1) calculated using eq. 
(13) are shown by the thick dotted line. The open square is 
lim n ^oo A c (n, n — 1; 5 C ). 




Fig. 6. Ratio r n = A„/A n -i plotted against 1/n. The solid line 
is the least-squares fit by the second-order polynominal of 1/n. 



A 8 = 116.9402383, 



(15) 



up to n = 8. The ratio r„ = A n /A n _i is plotted against 
1/n in Fig. 6. This ratio tends to converge to ~ 1.76 
and does not oscillate with n. Hence, we speculate that 
A n > for all n. This implies that an infinite series of 
phase transitions take place for 5 ~ 1. 

5.2 S = 5 C 

Because the infinite series of phase transitions are most 
prominent around S = S c in Fig. 5, we investigate this 
point in more detail. At this point, the ground state of 
the infinite AAFH1 is on the Gaussian critical point de- 
scribed by the conformal field theory with a conformal 
charge c =1. The ground state of a cluster- n is also de- 
scribed as a finite-size AAFH1 at the Gaussian critical 
point. It is known that the finite-size ground-state energy 



of the Gaussian model with an open boundary condition 
behaves as 29 ) 

E G (L,S C ) = Li G (oo,S c ) + e h + ^, A=^, (16) 

for large L, where L is the length of the chain, is the 
boundary energy, v is the spin wave velocity, and c (= 1) 
is the conformal charge. 

Fitting the ground-state energy obtained using the 
finite-size DMRG for 21 < L < 55, we obtain 



? G (oo, S c ) = -1.429085 ± 0.000005, 
e b = 1.1045 ±0.0005, 
A = 0.41 ±0.005, 



(17) 



where errors are estimated by changing the smallest sys- 
tem size used for the fitting from 21 to 29. 
Substituting eq. (16) in cq. (7), we obtain 



A c (n - l,n;6 c ) 



-e G (oo,5 c ) +e b + A 
( An + 



An + 1 



4n z 



1 



~ 2.5335 + 0.41 



\An 2 - 1 

The width of the DCn phase is given by 

8(n + l) 



(18) 



AX(n;5) =A 



(2n- l)(2n+l)(2n + 3) 
8(n+l) 



0.411- ' ^- . (19) 

V(2n- l)(2n+ l)(2n + 3)/ 1 ' 

This expression is positive for all n > 1. Therefore, an 
infinite series of phase transitions take place at A c (n — 
l,n;6 c ). It should be remarked that the width of the 
DCn phase decreases with n algebraically, while it de- 
creases exponentially for 5 ~ 1 according to eq. (14). 
This explains why the infinite series of transitions are 
most visible at S = S c in Fig. 5. 

We can also estimate A c (n — 1, oo; 8 C ) as 

A c (n - 1, oo; <5 C ) = -e G (oo, S c ) + e b + — - ^ , (20) 

by substituting eq. (16) in eq. (11). Thus, we obtain 

2nA 



X c (n - 1, n; S c ) - A c (n - 1, co; S c ) 



(21) 



4n 2 - 1 ' 

This implies that A c (n — l,n;5 c ) > A c (n — l,oo;<5 c ) for 
all n > 1. This confirms our numerical conclusion that 
no direct transition from the DCn phase to the DCoo 
phase takes place. The critical values A c (n — 1, n; S c ) arc 
plotted against n in Fig. 7. 

5.3 Physical quantities 

Physical quantities in the ground state is estimated 
in the same way as in the case of the undistorted 
MDC. 20 ' 21 ) Because each cluster-n carries a spin 1 degree 
of freedom, the low temperature magnetic susceptibility 
obeys the Curie law with the Curie constant C per unit 
cell given by 

°-WPry < 22 > 
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Fig. 7. Infinite series of critical points A c (n — l,n;<5 c ) at S = 
S c . Open circles arc calculated with eq. (7) using the values of 
Ea(2n + 1,5 C ) obtained by numerical diagonalization and finit- 
sizc DMRG. Filled circles are obtained using eq. (18). The open 
square is linin-s-oo A c (n, n — 1; <5 C ). The inset is an enlarged figure 
for large n. 




Fig. 8. A-dependence of Curie constant C at S = <5 C 

and the residual entropy per unit cell is given by 

1 In 3, 



n + 1 



(23) 



in the DCn phase. 

On the phase boundary between the DC(n — 1) and 
DCn phases, the residual entropy is larger than those 
in the DCn and DC(n — 1) phases owing to the mixing 
entropy of cluster-n's and cluster-(n— l)'s. It is given by 

— , (24) 



S = - In 
where x is the solution of 



(n + 1) In x — n ln(l — x) = In 3. 



(25) 



Although these formulae (22), (23), and (24) are com- 
mon to all values of S and n including the uniform case 
(S = O), 20,21 ' the actual A-dependences of these physical 
quantities show an infinite series of steps, if an infinite 
series of phase transition take place. These are plotted 




Fig. 9. A-dependence of residual entropy S at S = 5 C . Open circles 
are values on the phase boundary. 



against A in Figs. 8 and 9 for S = S c , where such behav- 
ior is most clearly observed. These structures are remi- 
niscent of the devil's staircase observed in various frus- 
trated systems and quasiperiodic systems. 30 ) In the case 
of the devil's staircase, however, physical quantities are 
quantized to arbitrary rational numbers and the width 
of the plateau decreases as the denominator increases. In 
contrast, in present model, the quantized values of the 
Curie constant and residual entropy are proportional to 
l/(n + 1) and the numerator is fixed. 

6. Summary and Discussion 

The ground-state phases of the alternating bond MDC 
with spins 1 and 1/2 are investigated. Owing to the local 
conservation laws, the ground states are rigorously con- 
structed, once the ground states of AAFHls are known. 
Each ground-state phase is described as a DCn phase 
that consists of an uniform array of cluster-n's. A cluster- 
n is equivalent to the AAFH1 with a length of 2n + 1. 

The ground-state phase diagram in the parameter 
space of the frustration A and the distortion S is numeri- 
cally determined. For S — 5 C ~ 0.2598, it is also analyzed 
using the conformal field theory. For 5 ~ 1, the cluster 
expansion analysis is carried out. Because both analy- 
ses predict an infinite series of phase transitions with 
respect to A, wc speculate that they should take place 
in the entire range of S c < S < 1. The behavior of the 
Curie constant and residual entropy is shown to exhibit 
an infinite series of steps. For small A, the ground state is 
the DCoo state, which corresponds to the Haldane state 
for S < 5 C and to the dimer state for 5 > S c . 

Thus, we found that the introduction of an alternat- 
ing bond distortion gives rise to a rich variety of quan- 
tum phases and phase transitions in the MDC. Other 
types of distortion that violate the local conservation 
laws produce even a richer variety of phases. These will 
be reported in a separate paper ,where the physics of the 
predicted phenomena are very different from the present 
one. 

We also stress that the possibility of the experimental 
realization of the MDC substantially is increased by al- 
lowing the lattice distortion. We hope the experimental 
synthesis of the MDC and the observation of the exotic 
phenomena predicted in the present paper in the near 
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future. Variations in S and A would be induced by the 
application of pressure. The infinite series of quantum 
phase transitions would manifest themselves as the pres- 
sure dependences of the low-temperature susceptibility 
and residual entropy. The substitution of nonmagnetic 
constituent atoms would result in similar effects. Because 
of the massive degeneracy in the DCn ground states, the 
magnetically ordered states would appear in the presence 
of an interchain interaction. Each ordered phase should 
have the periodicity of the corresponding DCn phase in 
the chain direction. Such a magnetic structure is observ- 
able by neutron scattering experiment. 

The numerical diagonalization program is based on the 
package TITPACK ver.2 coded by H. Nishimori. The nu- 
merical computation in this work has been carried out 
using the facilities of the Supercomputer Center, Insti- 
tute for Solid State Physics, University of Tokyo and the 
Supcrcomputing Division, Information Technology Cen- 
ter, University of Tokyo. KH is supported by a Grant- 
in- Aid for Scientific Research on Priority Areas, "Novel 
States of Matter Induced by Frustration" from the Min- 
istry of Education, Science, Sports and Culture of Japan 
and by a Grant-in-Aid for Scientific Research (C) from 
the Japan Society for the Promotion of Science. KT and 
HS are supported by a Fund for Project Research from 
Toyota Technological Institute. 
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